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of the flow similarity assumption adopted in the vicinity of the wall and in the core/outer region 


Keywords: of the flow field. Hence, the wall skin friction could possibly be doubtful when predicted using the 
Turbulent flow Prandtl-von Karman relation, in particular at high Reynolds numbers compared with values obtained 
Friction law from direct measurements. We therefore briefly review the recent advances in the logarithmic friction 
Pipe relation based on the latest pipe, channel, and boundary layer wall friction data over a wide range of 
Channel Reynolds numbers. Common similarities and differences among those flows, particularly in terms of 


Flat plate boundary layer the mean flow and the wall skin friction data are highlighted. Recent pipe friction relations (Mckeon 


et al., 2005) are compared with the Prandtl-von Karman logarithmic friction law and verified with 
new accurate pipe friction data. A revisited logarithmic skin friction relation for plane-channel flows 
(Zanoun et al., 2009) is reported, predicting the wall friction with an accuracy of better than +1.44%. 
A modified logarithmic skin friction relation for external flows introduced by Nagib et al. (2007) is 
presented and discussed. We believe that having a concise review and summary of the often used 
friction relations in wall-bounded shear flows in a single article will be of interest to fluid mechanics 
readers. 
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1. Introduction turbulence studies. In engineering, long gas and oil pipelines mo- 
tivate the examination of the skin friction relations available for 

The wall skin friction is of vital consideration in fundamental 
research and in practical applications. In fundamental research, 
the friction velocity is an important scaling parameter in wall 


flows along hydraulically smooth and rough pipe wall surfaces. 
Moreover, the wall skin friction accounts for approximately 50% 
of the total drag of a long-range subsonic airplane, and for 30% for 
* Coiesponding aithor: ground vehicles [1]. Thus, accurate and preferably independent 
E-mail address: el-sayed.zanoun@b-tu.de (E.-S. Zanoun). measurements of the wall skin friction with high enough spatial 
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resolution are needed. A common method to achieve this in two- 
dimensional boundary layer flows is to relate the wall shear stress 
(Ty) to the strain rate (dU /dy): 

dU 
Ty = a 
where y is the wall-normal distance, U is the local mean velocity, 
and yw is the fluid dynamic viscosity. This requires reliable ve- 
locity data within the viscous sublayer that are free from wall 
effects [2,3]. In most aerodynamics applications, however, the 
viscous sublayer is too thin to be resolved with high enough 
spatial resolution, making it difficult to obtain accurate velocity 
data within it, in particular, at high Reynolds number [4,5]. In 
spite of that, accurate measurements of the wall-shear stress 
were carried out, using the near-wall velocity data obtained by 
utilizing a micro-particle image velocimetry system with spatial 
resolution down to single-pixel ensemble correlation [4,5]. Here, 
one would conclude that if the size of the control volume of the 
measuring technique is not small enough to resolve the laminar 
sublayer, the so-called inertial sublayer might be used to estimate 
the wall skin friction [6]. Clauser [6] proposed a technique based 
on a claim that close to the wall, regardless of the flow geometry, 
a universal velocity function exists, ie. U = f(y, ty) [7], and 
extends outside the viscous sublayer, i.e. to the inertial sublayer 
and behaving in a logarithmic manner; more details are given in 
Section 3.3. 

In this paper, the logarithmic friction law is, therefore, dis- 
cussed for two classes of wall-bounded shear flows, i.e ducted 
flow (circular pipe and rectangular channel) and flat plate bound- 
ary layer of zero pressure gradient. They are fundamentally differ- 
ent flows in the sense that the pipe/channel flows are fully devel- 
oped whereas flat plate boundary layer are developing flows [8,9]. 
It is worth noting here that the logarithmic friction law in both 
flows has been determined based on a logarithmic scaling of the 
mean velocity profile at high enough Reynolds number. How- 
ever, it might incorrectly predict the skin friction if inappropriate 
values for the logarithmic law constants, i.e. the von Karman 
constant («) and the additive constant (B), are used [10,11]. A 
review of the logarithmic law constants showed variations of « 
by as much as 10%, which is considered a major source of error 
in determining the shear stress using the logarithmic friction law 
if « is not properly chosen [12]. Hence, to assess the logarithmic 
friction relations, new experimental data from the Cottbus large 
pipe, the CICLoPE pipe, and plane-channel facilities have been 
used. 

The importance of the logarithmic law of the wall in obtaining 
wall friction data motivated the present authors to highlight the 
mean velocity profile in Section 2. A summary of some available 
logarithmic friction relations for fully developed turbulent pipe 
flow is given in sub-Section 3.1. In sub-Section 3.2, the logarith- 
mic friction relation for plane-channel flow introduced by Zanoun 
et al. [13] is discussed and validated using various data sets from 
channel flows. New flat plate boundary layer (FPBL) wall skin 
friction data using the small wind tunnel facility at Brandenburg 
University of Technology (BTU-C-S) are presented and discussed 
in sub-Section 3.3. In addition, to extend the Reynolds number 
range for FPBL flows, the wall skin friction data from experiments 
at Illinois Institute of Technology (IIT) and Royal Institute of 
Technology (KTH) were used and are discussed. Conclusions and 
final remarks are given in Section 4. 


, (1) 


2. The mean velocity profile 


The law of the wall, U* = f(y*), was approached by Prandtl 
[7] using the dimensional analysis, while the velocity-defect law, 
U} — Ut = g(ņ), was proposed by von Kármán [14], where 
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Table 1 
Experimental values for the logarithmic velocity profile constants (« and B), 
equation Eq. (2), and the wake parameter (M). 


Flow K B IT Reference 
FPBL 0.384 4.127 0.50 [17] 
Pipe 0.39 4.42 0.23 [18] 
Channel 0.37 3.70 0.06 [12] 


Ut stands for either the normalized centerline pipe/channel local 
mean velocity or the freestream velocity in the case of flat plate 
boundary layer flow, U* is the normalized local mean velocity, 
and n = y/6, ô being the boundary layer thickness. On the other 
hand, in the inertial region, the logarithmic velocity profile was 
first introduced to the fluid mechanics community by von Kar- 
man [14] from similarity arguments, while using asymptotic anal- 
ysis by Milikan [15], and recently was derived by Oberlack [16] 
using the method of Lie-group symmetry analysis: 
Ut = ~Iny* +8, (2) 
where «x is the von Karman constant, and B is the additive con- 
stant, see Table 1. 

In a more general form, Coles [19] suggested the so-called law 
of the wake to describe the mean velocity distribution in both the 
overlap and outer regions: 


(3) 


assuming that the wake function w(7) is universal for all two- 
dimensional turbulent boundary layer flows which by definition 
w(0) 1 and w(1) 2, and M is called a profile or wake 
parameter. Fig. 1(a) presents the inner scaling of the mean ve- 
locity profile in two ducted flows, i.e. pipe and channel, and in a 
flat plate boundary layer for almost the same Reynolds number. 
The comparison in Fig. 1(a) was, therefore, made to examine 
the differences in the mean velocity profiles between the two 
ducted flows [9], versus the flat plate boundary layer [20] with 
reference to the logarithmic velocity profile. Plausible agreement 
between the types of flows is observed along the overlap region, 
however, this is not the case when the comparison is made either 
close to the wall or to the core/outer region. Thus Fig. 1(b) is 
produced, indicating that the experimental velocity data deviate 
from the logarithmic velocity profile for wall-normal distances 
y* < 150 and yt > 0.15R* where the velocity either overshoots 
or undershoots the logarithmic line. Values for constants of the 
logarithmic line summarized in Table 1 are used with respect 
to flow type to estimate the deviation of the measured velocity 
data from the logarithmic line and are presented in Fig. 1(b). 
Small differences, however, observable between the two types of 
flows within the viscous sublayer are obtained. On the contrary, 
the figure indicates discernible velocity differences between both 
types of flows in the core/outer region. Fig. 1(b) illustrates that 
the deviation from the logarithmic law in the core region of 
the channel flow is considerably less than those for the pipe 
and the flat plate boundary layer flows. To be more specific, 
the wake parameter, IT, which measures the departure of the 
mean velocity from the logarithmic velocity profile in the core 
region in both ducted flows and flat boundary layer flow was 
calculated. It is estimated to be 0.5 for the boundary layer flow 
case [20] presented and found to be in close agreement with the 
value obtained by Nagib et al. [17]. On the other hand, /7 for the 
pipe data was found to be 0.23, which is also in alignment with 
values obtained by Nagib et al. [17]. In a similar manner, the wake 
parameter in plane-channel flow was estimated and found to be 
0.06, which is much smaller than values obtained for the pipe and 
the boundary layer flows. Thus, one might conclude based on the 


1 II 
Ut = Inyt + B w(n), 
K K 
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Fig. 1. (a) The inner-scaled mean velocity distribution in pipe, plane channel, 
and flat plate boundary layer flows compared to the logarithmic velocity profile, 
and (b) the normalized velocity fractional deviation from the logarithmic line in 
correspondence with « and B values given in Table 1. 


relative magnitude of JT in those two types of the flows that the 
largest deviation was observed in the flat plate boundary layer, 
then the pipe flow, while the channel flow showed the smallest 
deviation. This might be attributed to a claim made by [17] that 
the presence of a favorable pressure gradient in duct flows tends 
to make JT smaller when compared with flat plate boundary 
layer flows. Hence, the logarithmic law may be valid much closer 
to the line of symmetry of the flow in the case of the channel 
flow and thus can be utilized to estimate the bulk flow velocity, 
resulting in reliable use of Eq. (6) and/or the logarithmic friction 
relation, i.e. equation (10), as will be shown in Section 3.2. On the 
other hand, the large deviation in the core/outer region from the 
logarithmic velocity profile in both the pipe and boundary layer 
flows might result in inaccurate estimation of the wall friction 
using the logarithmic friction relations. 
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3. The logarithmic friction law 
3.1. Pipe flow 


The logarithmic law of the wall provides a means for esti- 
mating the wall friction. Hence it is worth briefly summarizing 
the recent advances in the logarithmic friction relations in wall- 
bounded shear flows, starting with the pipe friction law. In the 
pipe shear flows, the general form of the logarithmic friction law 
was derived by Prandtl [7] in the form: 


/8/d = (1/«) In(Rep/A/32) + B — 3/2k, (4) 


based on a complete similarity assumption of the mean velocity 
profile in both the inner and the outer/core flow regions, where A 
is called the friction factor and Re» is the bulk Reynolds number 
defined based on the bulk flow velocity (Up), the pipe diameter 
(D), and the fluid kinematic viscosity (v). Prandtl [7] adopted a 
value of 0.4 for the von Karman constant («), and 5.5 for the 
additive constant (B) based on the analysis of Nikuradse’s mean 
velocity profiles in smooth pipe flow [21], resulting in the fol- 
lowing known Prandtl-von Karman logarithmic friction relation: 


1 
— =2 log(Re,v/7) — 0.8. 


Alternatively, based on a direct measure of the bulk flow ve- 
locity (U,) and the wall shear velocity (u+), the following friction 
factor (A) is commonly in direct use: 


ur 2 
nag =8( 3) , 
b 


where C; is the skin friction coefficient, u, is wall friction velocity 
defined as u, = ./Ty/p, p is the fluid density, and t, is the wall- 
shear stress that can be obtained either by using mean pressure 
gradient measurements or directly utilizing, for instance, oil film 
interferometry (OFI) [12,20,22]. Hence, it is worth noting here 
that if U, is not precisely obtained it is a source of significant error 
in estimating the skin friction data. 

In 1913, Blasius [23] introduced the one-seventh power law 
variation of the mean velocity to the fluid mechanics community, 
and used it to obtain the following power skin friction relation: 


(7) 


The upper limit of the Reynolds number proposed by Blasius [23] 
for using Eq. (7) is Re, = 10°. It might be considered, however, as 
the lower limit for the applicability of Eq. (5) [11,24]. One might 
also note that the power relation is in-line with the proposed 
Reynolds number-dependent power law representation of the 
mean velocity profile by [15,25,26]. 

On the other hand, based on the recent pipe friction data from 
the Superpipe facility, McKeon et al. [10] proposed the following 
pipe friction relation: 


1 
fe = 1.930 log( Rev) — 0.537, 


predicting the pipe friction for the Reynolds number range 31 x 
104 < Re, < 18 x 10°. 

Utilizing the new pipe facilities at BTU Cottbus-Senftenberg, 
and Bologna University, precise measurements of the bulk flow 
velocity (Up), and the mean pressure gradient (dP/dx) along the 
pipe test sections were carried out in CoSmaPipe [27,28], Co- 
LaPipe [29] and CICLoPE [30], resulting in the pipe wall friction 
data presented in Fig. 2, independently of the mean velocity 
profile measurements. Fig. 2 depicts a comparison among the 


(5) 


(6) 
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Fig. 2. The pipe friction factor, A, versus the bulk Reynolds number, Re), from 
pipe facilities in fully developed turbulent flow regime, (a) 3 - 104 < Re, < 10°, 
(b) 3-10 < Rep < 2- 10°. 


various data sets of the pipe friction versus the pipe friction rela- 
tions discussed in the present paper. An overall view shows good 
agreement among all experimental data sets, however, further 
discussion is warranted. Fig. 2(a) illustrates satisfactory agree- 
ment of the present experimental data for Rep < 10° with the 
Blasius friction relation, i.e. equation (7), within an accuracy of 
+1.08% with reference to Eq. (7), excluding the first three data 
points that deviate by more than 5%. Fig. 2(b) presents data for 
Re, > 10° compared with the various logarithmic friction rela- 
tions discussed in this paper. Eq. (5) predicts the CoLaPipe wall 
friction data with a maximum deviation of +1.124%, while the 
maximum deviation becomes +1.56% if equation (8) is consid- 
ered. On the other hand, the pipe friction data from the CICLOPE 
facility deviate by +1.537% and +0.51% with respect to Eqs. (5) 
and (8), respectively. 

Further, to correct for the deviation of the experimental data 
from the logarithmic friction law in the pipe flow due to the vis- 
cous effects, McKeon et al. [10] introduced the following friction 
relation: 


ie = 1.920 log( Res) — 0.475 — 7.04/(Rep, V/A)". (9) 
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Fig. 3. Representation of the logarithmic and the power friction relations 
discussed in this paper. 


Eq. (9) better predicts the pipe friction for the low range of the 
Reynolds number Re, < 10° since it accounts for the wall viscous 
effects. It was concluded by McKeon [10] that Eq. (9) predicts the 
friction factor to within 1.4% of the Princeton data (0.6% at high 
Reynolds number, 310 x 10? < Re, < 30 x 10°) and within 2.0% 
of the Blasius relation at low Reynolds numbers (10 x 10? < 
Re, < 90 x 10%). The data sets from Zanoun et al. [28] agree 
within +0.69% when compared with Eq. (9). On the other hand, 
the CoLaPipe friction data for 16 x 104 < Re, < 10° agree within 
1.65% compared to Eq. (9), while the CICLoPE pipe wall friction 
data showed better agreement, i.e. within +0.5%, with Eq. (9) for 
13 x 104 < Re, < 2 x 10° (see Table 2). 

To evaluate all friction relations discussed earlier, Fig. 3 repre- 
sents a comparison of the logarithmic friction relations, i.e. Eqs. 
(5), (8) and (9), and the power friction law, Eq. (7). From the 
figure, one hardly observes any differences among all the rela- 
tions, however, the deviations of McKeon et al.’s [10] relations 
(i.e. Eqs. (8) and (9)) with respect to the Prandtl-von Karman 
logarithmic law, i.e equation (5), and Blasius power relation (7) 
are discernible, as Fig. 4 illustrates. The difference of both log- 
arithmic relations, i.e. Eqs. (8) and (9), introduced by McKeon 
et al. [10], from the Prandtl-von Kármán law increases as the 
Reynolds number increases. Within the present working range 
of the Reynolds number, i.e. 6 x 104 < Re, < 2 x 10°, of both 
the CoLaPipe and the CICLoPE facilities, the deviation of both 
Eqs. (8) and (9) from the Prandtl-von Karman law is approxi- 
mately 2%, however, it increases with increasing Reynolds num- 
ber, for instance, it reaches 3% at Rey = 2 x 10’. Moreover, 
the friction relations of McKeon et al. [10] deviate within +2% 
for 104 < Re, < 10° with respect to the Blasius power friction 
relationship, however, Eq. (9) behaves better for the low range 
of the Reynolds number. 

Table 2 briefly summarises the calculated deviations in the 
present pipe friction measurements compared to the pipe friction 
relations discussed in text. 


3.2. Plane-channel flow 


In plane-channel flows, the two-dimensionality of the flow 
is expected for channels having high enough aspect ratios, i.e. 
width(b)-to-full height(2h) ratio AR = b/2h > 7 [31]. In addition 
to flow two-dimensionality, the fully developed turbulent state in 
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Table 2 
Summary of calculated deviations in the present pipe friction measurements compared to the pipe friction relations discussed in 
text. 
Facility Equation (7)[23] Equation (5)[7] Equation (8)[10] Equation (9)[10] 
Re, < 10° Re, > 10° 31-104 < Re, < 18-107 31-104 < Re, < 6-107 
10-10? < Re, < 9-104 
CoSmaPipe +1.081% = +0.50% +0.69% 
CoLaPipe +1.424% +1.124% +1.56% +1.65% 
CICLoPE = +1.537% +0.51% +0.5% 
9.0x10° T T T T T T T 
Deviation of equation (8) from equation (5) === Channel (Ch) log. relation: equation (10) [13] | | 
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—-— Deviation of equation (8) from equation (7) 8.0x10° -+ @ Coil film interfeormetry, Ch b/2h=12 [13] H 
— — - Deviation of equation (9) from equation (7) z C, pressure gradient, Ch b/2h=11.7 [33] 
© Channel, b/2h=10 [34] 
R 3 A Channel, b/2h=18 [2] 
7.0x10" = --e-- Pipe, log. relation: equation (6) [10] 
x Cc —-—Channel log. relation, x=0.41, B=5.17 [31] 
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Fig. 4. Deviation of McKeon et al. [10] friction relations from the Prandtl-von 
Karman [7] and the Blasius [23] friction relations. 


channel flows is assured for long enough entrance lengths [12]. 
Thus, the local wall skin friction, Cf = 2(u,/Up)’, data in plane- 
channel flow was obtained using either the mean pressure gra- 
dient or oil film interferometry where flow was assured to be 
fully developed. The evaluation of C; was then established based 
on both measuring techniques and the results obtained are pre- 
sented in Fig. 5 in the form of C as a function of Re,, where Rep 
is the Reynolds number defined based on the bulk flow velocity 
(Up), the channel full height (2h), and the fluid kinematic viscosity 
(v). The wall skin friction data obtained using a channel that 
had an aspect ratio of 12, from the mean-pressure gradient, was 
observed to compare well with the local skin friction measured 
by oil film interferometry [12]. Thus, based on those experimental 
skin friction data in a plane channel for 17 x 10? < Re, < 
25 x 10* (900 < R* < 5053), Zanoun et al. [13] proposed the 
following revisited logarithmic skin friction relation: 


1 


Ja 
noting that it was obtained utilizing the logarithmic law of the 
wall by adopting a value of 0.37 for the von Karman constant («), 
and 3.7 for the additive constant (B). It was referred to, recently, 
as the Prandtl-Zanoun relation by June et al. [32] taking the 
following modified form: 


= 1.917 In(Re,/G) — 1.286, (10) 


= 2.2 logy(RepWA) — 1.9656 


Vi 
where A = 4C;. They [32] defined Rep as the Reynolds number 
based on the hydraulic diameter (Dy), i.e. Di = 4A/P, A is the 
cross-sectional area of the duct, and P is the wetted perimeter. 
The logarithmic friction relation, i.e. equation (10), is pre- 
sented in Fig. 5 as a function of the mean-based Reynolds number, 


(11) 
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Fig. 5. Representation of the logarithmic skin friction relation, Eq. (10), 
compared with recent experimental results and data extracted from the 
literature. 


Rep, compared with recent channel experimental results and 
channel data extracted from the literature. Error analysis using 
the least-squares method indicates that the channel experimen- 
tal skin friction data deviate within +1.44% when compared 
with Eq. (10). Such good agreement obtained might be attributed 
to the fact that the logarithmic mean velocity profile in the core 
region of the channel flow deviates much less than those for 
pipe flow, as well as for flat plate boundary layer flow as shown 
earlier in Fig. 1(b). This mean velocity behavior results in an 
accurate estimation of the bulk flow velocity, and hence, precise 
estimation of the wall skin friction data using either equation 
(6) and/or equation (10). Fig. 5 illustrates also the skin friction 
data obtained from mean pressure gradient measurements for 
4x 104 < Re, < 19x 104 (10? < Rt < 4x10?) ina plane channel, 
having an aspect ratio of 11.7:1 [33]. The following logarithmic 
friction relation was then proposed by Monty [33] based on his 
own channel experimental skin friction data: 


ma = 4.175 log,9(Rep./c) — 0.416, 
f 

which is also presented in Fig. 5. Eq. (12) agrees with Eq. (10) 
within an accuracy of +0.7% for Re, < 1.2 x 10°, however, its 
deviation increases with increasing Reynolds number, e.g. reach- 
ing 1.3% for Re, ~ 2.5 x 10°, see Fig. 6. On the other hand, a 
logarithmic friction relation employing the constants proposed by 
Dean [31], i.e. x = 0.41 and B = 5.17, exhibits less agreement 
with the channel experimental data and with both friction rela- 
tions, i.e. Eqs. (10) and (12), in particular, for the high range of the 
Reynolds number. The low Reynolds number data from [34] and 
from [2] are also presented in Fig. 5, showing satisfactory agree- 
ment with Eq. (10). For contrast, and to emphasize the difference 


(12) 
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between the fully developed turbulent channel and pipe flows, 
we have included the logarithmic friction relation, i.e. equation 
(8), for the pipe flow by McKeon et al. [10] which extends to 
very high Reynolds number, 310 x 10° < Re, < 18 x 10°. On the 
other hand, the deviation of both power friction relations pro- 
posed by Dean [31], G = 0.73 Re. and Zanoun et al. [13], 
C = 0.743 Re, ^”, from Eq. (10) are plotted in Fig. 6. One might 
note that the power skin friction relation proposed by Dean [31] 
agrees better with Eq. (10), however, for a narrow range of the 
Reynolds number (6.5 x 10? < Re, < 15 x 104) when compared 
with the power friction relation proposed by Zanoun et al. [13]. 


3.3. Boundary layer flow 


Along flat plate boundary layer flows, the wall skin friction (Cp) 
can be found by direct measurement of the wall friction velocity 
(u), using, e.g., oil film interferometry (OFI) [22]. Alternatively, it 
can be obtained using Eq. (1) and/or by calculating the gradient of 
the momentum thickness (@) of the boundary layer with respect 
to the streamwise distance (x), i.e. d0 /dx, based on the Karman 
integral momentum equation: 

Ur 


CG do z 
Us) ` 


2 da 

where U» is the free stream mean velocity. The estimation of 
dd /dx, however, requires measuring the streamwise mean veloc- 
ity versus the wall-normal distances, in particular, in the vicinity 
of the wall at various locations along the streamwise direction 
of the boundary layer. It seems at first glance that it is easy 
to measure the local mean velocity versus the wall-normal dis- 
tance, however, in most aerodynamics applications, the viscous 
sublayer is too thin to be resolved with high enough spatial 
resolution [4,5], making it difficult to obtain accurate velocity 
data within it, particularly, at high Reynolds number. Hence, if the 
control volume of the measuring technique is not small enough 
to spatially resolve the viscous sublayer, it is proposed that the 
logarithmic region of the mean velocity profile should be utilized 
to obtain the wall skin friction. Clauser [6], hence, introduced the 
following plot based on the logarithmic velocity profile: 


Cr yur Cr 
= Z alo (=) A logo ( # \ | +B, 
U 5 | Si0( >) + S10 ( z j| t 


where U is the local mean streamwise velocity component, A 
=ln(10)/x, « is the von Kármán constant, and B is the additive 
constant. It is, however, worth re-noting that the use of the 
Clauser plot is linked to the logarithmic velocity profile within 
the inertial sublayer in boundary layer flows [6]. To obtain a 
value for the wall skin friction utilizing the Clauser method, it 
is usual to superpose the Clauser plot, i.e. equation (14), on a 
measured velocity profile within the logarithmic range of the 
inertial sublayer via the adjustment of Cp as presented in Fig. 7. A 
proper value for the skin friction coefficient can then be obtained 
when the Clauser plot matches well the measured mean velocity 
profile within the logarithmic range, see Fig. 7. However, it is 
obvious from Eq. (14) that there is a direct dependence of the 
skin friction coefficient upon the values adopted for the von 
Karman constant (x ) and the additive constant (B). Crook [35], for 
instance, concluded that the prediction of u, is sensitive to the 
values of « and B adopted in the logarithmic law. He observed 
that a +0.5% change in the slope of the logarithmic line, 1/x, 
results in a 12% difference in u,. Hence, a careful choice for the 
slope of the logarithmic line is an important issue in obtaining 
the wall skin friction using the logarithmic friction law. When the 
most recent and accurate values for both constants A = 5.996 
(x = 0.384) and B = 4.127 obtained [17,20] from flat plate 


(13) 


a (14) 
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Fig. 6. Deviation of the experimental data and channel friction relations dis- 
cussed in the paper from the logarithmic skin friction relation, Eq. (10), proposed 
by Zanoun et al. [13]. 


boundary layer measurements are adopted, they would result in 
a rather good application of the Clauser method. It is observed, 
however, that the recent values for both « and B obtained [17,20] 
are in very close agreement with values « = 0.3806 (A = 6.05) 
and B = 4.05 obtained by careful measurements carried out by 
Winter and Gaudet [36]. On the other hand, the wall skin friction 
(Cy) results obtained utilizing equation (14) with recommended 
constants [17,20] were found to be in fair agreement with Cp data 
obtained when the most often used values A = 5.616 (x = 0.41) 
and B = 5.0, proposed by Coles and Hirst [37], for flat plate 
turbulent boundary layers were used. 

Fig. 8 presents the wall skin friction data from two measuring 
techniques, laser-Doppler anemometer (LDA) and oil film inter- 
ferometry (OFI), versus some friction relations extracted from the 
literature [17,38]. The figure illustrates satisfactory agreement be- 
tween the wall skin friction data obtained from the LDA near-wall 
mean velocity profiles and the OFI as well as with data obtained 
utilizing the Clauser plot. However, one might observe a little 
under-estimation of the wall skin friction using the OFI that might 
be attributed either to the accuracy of the measured oil viscosity 
or to the effect of the wall surface temperature at which mea- 
surements were carried out. A revisited F. Schultz-Grunow [39] 
logarithmic friction relation derived by Nagib et al. [17] is plotted 
in Fig. 8 and compared with the current experimental data. Good 
agreement with the modified F. Schultz-Grunow [39] relation as 
well as with the power friction relations [38] is obtained. One also 
observes from Fig. 8 the early transition to the turbulent regime at 
Re, ~ 10° by the artificial triggering of the boundary layer versus 
the natural transition at approximately Re, ~ 5 x 10°, where 
Re, is defined based on the streamwise freestream mean velocity 
(Ux), the streamwise distance (x), and the kinematic viscosity (v). 


Further examination, however, of the wall skin friction (Cp) 
along flat plates for high Reynolds numbers has been carried 
out, recently, confirming the necessity for direct and independent 
measurement of the wall shear stress (tẹ) [17]. Nagib et al. [17] 
found that many of the available empirical friction relations 
might describe well the local Cp behavior when modified and 
supported by accurate experimental data. Thus, they used the 
direct wall shear stress measurements of the experiments at 
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Fig. 7. The Clauser plot superposed on a selected case from the mean velocity 
experimental results obtained using a laser-Doppler anemometer. 
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Fig. 8. Recent wall skin friction data compared with modified laminar and 
turbulent formulae extracted from [17,38]. 


Illinois Institute of Technology (IIT), Royal Institute of Technology 
in Stockholm (KTH) and the National Diagnostic Facility (NDF) to 
study the development of Cp and evaluate the quality of the var- 
ious friction relations found in the literature. It was commented 
by Nagib et al. [17] that none of the original friction relations 
presented in their paper [17] were in good agreement with the 
recent experimental data obtained from KTH and NDF. They 
attributed the disagreement to the low Reynolds number effects 
to establish appropriate values for the coefficients required for 
the various forms of the skin friction relations. To revisit, for 
instance, a logarithmic friction relation for flat plate boundary 
layer flow, Nagib et al. [17] adopted Coles-Fernholz’s friction 
relation [40]: 


CG; = 2[1/« In(Reg) + BI? , (15) 


for further theoretical developments. Satisfying asymptotic con- 
sideration, Nagib et al. [17] proposed new numerical values for 
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Fig. 9. The boundary layer wall skin friction, Cp, versus the Reynolds number, 
Reg: 

(a) KTH and NDF experiments versus some selected friction relations [17], 

(b) deviation of KTH and NDF experimental data and Eqs. (16) and (17) 
from Eq. (15). 


both « = 0.384, and B = 4.127 based on the collective NDF and 
KTH oil-film data. The modified version of Coles-Fernholz friction 
relation is presented in Fig. 9a and further discussed. In addition 
to Eq. (15), Nagib et al. [17] modified F. Schulz-Grunow [39] 
original friction relation to read as: 


C; = 0.3475 [log(Rex) — 0.407764. (16) 


where Re, is defined based on the freestream mean velocity (Ux), 
the streamwise distance (x) measured from the leading edge, and 
the kinematic viscosity (v). Moreover, Nagib et al. [17] modified 
the empirical coefficients of the Prandtl-Kármán skin friction 
relation to be written as follows: 


G"? = 4 [log(Rexy C) — 2.12. 


Noting that the original and modified empirical coefficients in 
Eqs. (15)-(17) are documented in Table 1 in Nagib et al. [17]. To 
be able to compare both Eqs. (16) and (17) with Eq. (15) in Fig. 9a, 
Rex is related to Reg through the following revisited empirical 


(17) 
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relation: 


Reg = 0.012777 - Re? 8678 (18) 


which predicts Reg with an accuracy of +1% when compared 
with the experimentally determined Reg [17]. All three selected 
and revisited friction relations, i.e. Eqs. (15)-(17) are plotted in 
Fig. 9(a), showing good agreement within the experimental accu- 
racy. However, we find from Fig. 9(a) that the revisited version 
of the nearly a century old Prandtl-Karman skin friction relation 
appears to be exceptionally in good agreement with Eq. (15). 
The good collapse of skin friction relations corrected by Nagib 
et al. [17] and the KTH and NDF experimental data, independent 
of how the flow was triggered along the flat plates and how the 
Reynolds number was achieved presented in Fig. 9(a), indicates a 
self-similarity of the wall skin friction [17]. 

Albeit the wall skin friction data from both NDF and KTH were 
directly obtained and used to evaluate the quality of the various 
friction relations, their deviations from Eq. (15) were estimated 
and presented in Fig. 9(b). The deviation of the NDF data from 
the revisited Coles-Fernholz friction relation, i.e. equation (15), 
was found to be within +1.94%, whereas the KTH deviates within 
+ 1.29%. Fig. 9(b) indicates also that Eq. (16) deviates within 
+2.5% when compared with Eq. (15) for 1.2 x 104 < Rey < 10°, 
however, the deviation increases for 1.2 x 104 > Reg > 10°. 
On the other hand, the revisited version of the Prandtl-Karman 
friction relation appears to be in outstanding agreement with 
Eq. (15) over a wide range of Reg. 


4. Concluding remarks 


Different methods to determine the wall friction in turbulent 
pipe, channel, and flat boundary layer flows have been reviewed. 
In spite of the fact that some of the methods are still under dis- 
cussion the comprehensive presentation of the data may help to 
understand better the differences and limitations of the available 
logarithmic friction relations under various boundary conditions. 
A surprising outcome, however, of the present study is that the 
various data sets collected and illustrated here, fitted very well 
among each other and with the logarithmic friction relations 
presented. This, partially, might be due to the proper values 
chosen for constants of the logarithmic velocity profile, the high 
quality of the new experimental facilities and accuracy of the 
measuring techniques used. 

In closing, from the results and discussions, we conclude that 
precise and independent wall skin friction data are indeed of 
vital importance in engineering applications, and as a scaling 
parameter for the analysis of wall turbulence. 

The turbulent shear flows in closed conduits, i.e. circular pipe 
and rectangular channel, and boundary layers of zero pressure 
gradient are two classes of flows that are fundamentally different 
in the sense that pipe/channel flows are fully developed whereas 
boundary layers are developing flows; however, a comparison 
has been be made in terms of the mean velocity profiles and 
skin friction relations, reflecting the combined effects of flow 
geometry and pressure gradient. We conclude based on esti- 
mating the magnitude of the wake parameter, IT, in these two 
types of flows that the largest deviation in terms of the mean 
velocity in core region was observed in the flat plate boundary 
layer, then the pipe flow, while the channel showed the smallest 
deviation. This as claimed by Nagib et al. [17] might be attributed 
to the combined effects of flow geometry and favorable pressure 
gradient in duct flows. This observation in fact motivates further 
research in connection with clarifying the structural behavior of 
the outer layer between these two types of wall-bounded shear 
flows. 
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The Prandtl-von Karman logarithmic friction law with x = 0.4 
and B = 5.5 applies for the pipe flows with reasonable accuracy. 
The « value chosen is in-alignment with a recent remark made by 
Baumert [41] based on theoretical considerations that suggested 
the von Kármán constant « = 1//22 ~ 0.399 which is also 
very close to the «-value Nikuradse [21] suggested from his 
measurements in 1932. In agreement with McKeon et al. [10], 
the results suggest that the pipe friction factor behavior falls 
into three regimes. For Re, < 10°, the data are well repre- 
sented by the Blasius relationship, Eq. (7). For 10° < Re, < 
3.2 x 10°, Prandtl-von Karman’s friction law, Eq. (5), applies with 
reasonable accuracy. For Re, > 3.2 x 10°, the friction relations 
proposed by McKeon et al. [10] yield a better description of the 
pipe wall friction up to the highest Reynolds number, Re, ~ 
35 x 10°, achieved by Zagarola and Smits [24]. 

For plane channel flows, having a large enough aspect ra- 
tio, the logarithmic skin friction relation proposed by Zanoun 
et al. [13] with « = 0.37 and B = 3.7, up to more than twice 
the highest Reynolds numbers available, predicts the wall skin 
friction with high accuracy, in particular for high Reynolds num- 
bers. For contrast, and to emphasize the differences between the 
fully developed turbulent channel and pipe flows, the logarithmic 
friction relation, i.e. equation (8), for the pipe flow by McKeon 
et al. [10] which extends to very high Reynolds number was 
presented in Fig. 5, showing a clear difference between the two 
ducted flows in terms of the wall skin friction. 

A revisited logarithmic skin friction relation, i.e. equation (15), 
with properly chosen parameters, x 0.384 and B = 4.127, 
by Nagib et al. [17] is accurate in representing the wall skin 
friction in zero pressure gradient boundary layer experiments. 
In addition, the revisited version of the Prandtl-von Karman 
friction relation, i.e. equation (17), appears to be exceptionally 
in outstanding agreement when compared with Eq. (15) over 
a wide range of Rey. Finally, since the wall roughness strongly 
influences the above-discussed friction relations in this paper, it 
is of crucial importance to consider the roughness effect in an 
additional review. 

It is worth noting here that the research work, published in 
this paper, required collaborations and the first author is very 
happy that he had the opportunity to work together with C. 
Egbers, H. Nagib, F. Durst, G. Bellani and A. Talamelli to get 
reliable measurements in pipe, channel and flat plate boundary 
layer flows. The pipe flow data were obtained using the Co- 
LaPipe at BTU C-S by E.-S. Zanoun and C. Egbers, and CICLoPE 
at University of Bologna by G. Bellani and A. Talamelli. G. Bel- 
lani and A. Talamelli designed and prepared the CICLoPE set-up 
and E.S. Zanoun had an opportunity to work with both during 
one of the EuHit (European High-performance Infrastructures 
in Turbulence) campaigns in CICLoPE. The channel flow data 
were obtained using plane-channel facility at LSTM-Erlangen, 
Friedrich-Alexander-University, during the research stay of the 
first author working under the supervision of both F. Durst and 
H. Nagib. H. Nagib provided the flat plate boundary layer via his 
close collaboration with both KTH (Royal Institute of Technology 
in Stockholm) and NDF (National Diagnostic Facility). All authors 
contributed in writing, interpreting and discussing the data in the 
different sections of the manuscript. 
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